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Abstract 

The results on cosmic rays detected by the gravitational antenna NAUTILUS have 
motivated an experiment (RAP) based on a suspended cylindrical bar, which is 
made of the same aluminum alloy as NAUTILUS and is exposed to a high energy 
electron beam. Mechanical vibrations originate from the local thermal expansion 
caused by warming up due to the energy lost by particles crossing the material. The 
aim of the experiment is to measure the amplitude of the fundamental longitudinal 
vibration at different temperatures. We report on the results obtained down to a 
temperature of about 4 K, which agree at the level of ~10% with the predictions of 
the model describing the underlying physical process. 
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1 Introduction 



Cosmic rays generate signals in a massive gravitational wave (GW) detector. 
The signals are due to vibrations produced by the heating along the particle 
trajectories, depending on the thermal expansion coefficient and on the specific 
heat of the material. 

For cylindrical detectors the acoustic model [1,2,3,4,5,6,7,8,9] gives quanti- 
tative predictions for the amplitude of the longitudinal modes of oscillation 
due to the energy loss of a particle impinging on the cylinder. The ampli- 
tude is proportional to the ratio of the thermal expansion coefficient to the 
heat capacity. This ratio is part of the definition of the material Griineisen 
parameter: 



where (3 is the volume expansion coefficient, K T is the bulk module at constant 
temperature, V m is the molar volume and c v is the specific heat at constant 
volume. The parameter 7 is nearly constant over a wide range of temperatures. 

Cosmic ray measurements made with the aluminum GW detector NAUTILUS 
gave essentially two types of results: 

• When the NAUTILUS bar was at a temperature T=0.14 K, that is in a 
superconductive regime, the rate of the high energy signals due to cosmic 
ray showers was larger by two-four order of magnitude than the expectations 
based on the models [10,11] 

• When the bar temperature was T=1.5 K, aluminum in a normal state, few 
signals were large, most of the signals obeying [12] the predictions. 

It was not clear to us whether that depended on the behavior of the Griineisen 
parameter in the superconductive regime (for example, data on the aluminum 
critical field as a function of pressure and temperature [13] suggest values of 
the parameter larger than in the normal state [14]). 

To gain information on this problem a new experiment (RAP) was designed, 
consisting in bombarding a small aluminum bar, whose temperature we could 
control, with a known electron beam, provided by the Beam Test Facility 
(BTF) of the DA$NE ^-factory complex at the INFN Laboratory in Frascati. 

We compare the measured maximum amplitude of the bar fundamental lon- 
gitudinal mode of oscillation with the expected value: 
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where [3] : 



Bo = ~W (2) 
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Here a is the linear thermal expansion coefficient {(3 = 3a for aluminum), 
L, M are respectively length, mass of the cylinder and W is the total energy 
released by the beam to the bar. The term e accounts for corrections estimated 
by Monte Carlo methods due to contributions 0[(R/ L) 2 }, R being the cylinder 
radius, and to the beam structure (see Appendix A). 

After a short description of the experimental setup and a discussion on the 
beam simulated interactions within the bar, we report on the results obtained 
down to a temperature of ~ 4 K. 



2 Experimental setup 



The RAP setup is composed by the beam [15] and the detector [16,17]. The 
cylindrical test mass, the suspension system, the cryogenic system, the trans- 
ducer and the data acquisition system are parts of the detector. 



2. 1 Beam 



DA$NE BTF provides the controlled beam for RAP. Although the facility 
is optimized for the production of single electrons (or positrons), mainly for 
high energy detector calibration and testing purposes, it can provide particles 
in a wide range of energy (25-750 MeV) and intensity (up to ~ 10 10 parti- 
cles/pulse). Particles can be produced in pulses of 1 or 10 ns duration, with a 
fairly uniform distribution. The maximum energy of the beam is ~ 510 MeV 
when BTF is jointly operated with the injections into the DA$NE collider. 
The installation site of RAP is 2.5 m far away from the exit of the beam (see 
Fig. 1). The position of the beam at the exit of the line and at the entrance of 
the cryostat was monitored shot-by-shot by two high sensitivity fluorescence 
flags, 1 mm thick alumina doped with chromium targets. The spot size at the 
entrance of the cryostat, 50 cm far from the center of the bar, is ~ 2 cm in 
diameter. The Monte Carlo simulation of the detector (see Section 4) indicates 
that the beam spot at the surface of the cylindrical bar mantains almost this 
dimension, due to negligible effects induced by the thickness of the cryostat 
vacuum shields intercepted by the beam. 

BTF delivers to RAP single pulses containing O[10 8 ] electrons each. In this 
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Fig. 1. Schematic layout of the RAP setup in the BTF hall. 

range standard beam diagnostics devices to monitor the beam charge have 
low sensitivity, while saturation affects particle detectors. Therefore a moni- 
tor (BCM), based on an integrating current transformer readout by a charge 
digitizer, was installed to measure the pulse charge. The device accuracy is 
~3% and the sensitivity, dominated by the digitizer, is ~ 1.4 - 10 7 electrons, 
corresponding to a a ~ 4 • 10 6 electrons. The device is equipped with a cal- 
ibration coil used to control the gain, the noise and the time shaping of the 
generated signals. 



2.2 Test mass and suspension system 

The oscillating test mass is a cylindrical bar (diameter 0.182 m, length 0.5 m, 
mass 34.1 kg) made of A15056, the same aluminum alloy (5.2 Mg%, 0.1%Mn, 
0.1% C) used for NAUTILUS. The resonance frequency of the first longitudinal 
mode of vibration is 5096 Hz at T=296 K. 

The suspension is an axial-symmetric system of seven attenuation stages in 
series made of copper. The system provides a -150 dB attenuation of the 
external mechanical noise in a frequency window spanning from 1700 to 6000 
Hz. The symmetry axis of the suspension passes through the center of mass 
of the cylinder. The suspension is linked to the lateral surface of the cylinder 
by means of two brass screws at the mid-section of the cylinder. 
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2.3 Cryogenic system 



The cryogenic system is based on a commercial cylindrical cryostat (3.2 m 
height, 1.016 m diameter) and a 3 He /L He dilution refrigerator (base tempera- 
ture=100 mK, cooling power at 120 mK=l mW), to be installed in the near 
future. Cryogenic operations start by the filling of the cryostat with liquid 
nitrogen, both in the nitrogen and in the helium dewars, followed by the 
cool-down to 4.2 K with liquid helium. The thermal exchange among liquid 
helium and the bar is realized by filling the space surrounding the bar with ~ 
1 mbar of gaseous helium. This is indeed the only thermal link between the 
bar and the liquid helium, the suspension-bar assembly being mechanically 
disconnected from the cryostat, with the exception of three thin stainless steel 
wires connecting the two components. 

Temperatures are measured by two platinum resistor and three ruthenium 
dioxide resistor thermometers, in the 60-300 K range and in the 0-60 K range 
respectively, with a precision of 0.1 K. All the thermometers are controlled by 
a multi-channel resistance bridge. 

A mechanical structure encloses the cryostat allowing an easy positioning of 
the detector on the beam line and the consequent removal after the expiration 
of the dedicated periods of data taking. 

Fig. 2 shows the bar suspended inside the upper half of the cryostat. 




Fig. 2. The suspended bar in the opened cryostat. 
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2.4 Transducer and data acquisition system 

Two Pz24 piezoelectric ceramics (0.01x0.01x0.005 m 3 ), electrically connected 
in parallel, are inserted in a slot cut in the position opposite to the bar suspen- 
sion point and are squeezed when the bar contracts. In this configuration the 
stress measured at the bar center is related to the displacement of the bar faces. 
The generated voltage signal is sent to an amplifier (V no i se ~4 nV/vTIz). The 
data acquisition system (DAQ), based on a peak sensing 16-bit VME ADC and 
a VME Pentium III CPU running Linux, collects and stores on disk amplified 
data coming from the transducer and data originated from the thermometers 
and from the beam monitoring system. The acquisition system provides also 
tools for performing a fast inline analysis of the collected data. 



3 Calibrations 

It is possible to measure the value of A, the transducer electro-mechanical 
coupling factor, making use of the electric equivalent circuit shown in Fig. 3. 



This circuit includes &R\ — L\ — C\ series that expresses the bar characteristics 
(dissipation, mass, elasticity) and in parallel the Pz24's with capacity Ci- To 
measure the equivalent circuit parameters we proceed as follows: 

• We apply a signal V g cos(u t) at the Pz24's for a time At much smaller 
than the decay time of the oscillations, where u = 2irf = 1/ \JL\C\ is the 
angular resonance frequency of the bar. 

• We switch the Pz24's from the oscillator producing V g to the amplifier and 
measure the signal V Q cos(u t + </>) due to the bar oscillation. 




o 



Fig. 3. Equivalent electrical circuit and calibration scheme. 
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1.32 


71 


5397.3 
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1.48 


4.5 


5412.7 


84 


1.32 



Table 1 

Frequency, decay time of the first mode of longitudinal oscillation and transducer 
coupling factor as a function of temperature. 

It is possible to show [18] that C\ is given by the following equation: 
V 



VgAtirfo 



[C 2 + C L ] 



where Cx is the load capacity that includes the cable and the input impedance 
of the amplifier. Having measured the value of C\ we obtain A for the first 
mode of longitudinal oscillation by: 



A = 



2tt/ MCi 
C 2 + C L V 2 



(3) 



where M is the mass of the bar. 

The values of A, obtained at different temperatures and before the collection of 
data useful for the estimation of the vibrational effects induced by the beam, 
are shown in Table 1 together with the measured frequency of the first mode 
of longitudinal oscillation f a and the its decay time r Q . 

The PZ24's are inserted at a distance of ~ 0.09 m from the cylinder axis, 
while the calibration procedure we have used is correct in the (R/L) 2 <C 1 
approximation. We checked the calibration procedure inserting a commercial 
calibrated accelerometer with 5% accuracy at the center of one of the bar end- 
surfaces. At room temperature the displacement measured by the PZ24's was 
2.4% smaller than the one measured by the accelerometer and at T ~ 77 K 
the displacement was ~6% smaller. Consequently we assume a 6% systematic 
error in measuring the displacement amplitudes. 



4 Simulation of the beam energy loss in the aluminum bar 



We have developed a Monte Carlo simulation using the CERN package GEANT 
3.21, taking into account the real geometry and material of the cryostat and 
of the bar, and a realistic parametrization of the BTF beam, i.e. the beam 
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spot size and divergence, and simulating the passage of electrons through the 
matter. The size of the beam spot and its angular and momentum spread 
have been simulated according to the measured characteristics: the well fo- 
cussed BTF beam has a fairly Gaussian shape in the transverse plane with 
<J X = (T y = 0.005 m, the momentum spread is ~ 1% and the beam emittance 
is e = 1 mm mrad. The cryostat has been schematized as a succession of three 
coaxial aluminum cylinders. Due to the presence of these aluminum shields, 
with a total Al depth of 1.7 ■ 10~ 2 m, the energy of the electrons is slightly 
degraded by the emission of Bremsstralhung photons. On average, more than 
100 secondary particles are produced for each 510 MeV primary electron en- 
tering the bar. The spatial distribution of these secondary particles inside the 
bar is shown in Fig. 4 for the three perpendicular projections. The develop- 
ment of the electromagnetic shower is clearly visible; indeed the bar diameter 
corresponds to ~ 2 radiation lengths. 




-8 -4 4 8 

Z (cm) 



Fig. 4. Monte Carlo simulation: secondary particles distribution for 510 MeV pri- 
mary electrons impinging on the bar for the three projections. The black arrow in- 
dicates the direction of the incoming beam. 

The total energy released to the bar by an electron in the beam is then evalu- 
ated adding the contribution of the energy lost by all the JV, secondary particles 
(j = is the contribution of the primary electron itself), at each of the Mj 
tracking steps in the simulation code: 

j=0 k=l 

The distribution of the energy lost by one electron with energy of 510 MeV in 
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10 6 trials, as shown in Fig. 5, leads to an average energy loss per electron: 
< AE > ± a AE = 195.2 ± 70.6 MeV (4) 



and for a beam composed by N electrons to a total energy loss: 

W = N < AE > , a w = y/N a AE (5) 

We performed extensive simulations (see Appendix A) in order to: (a) apply 
the acoustic model to each infinitesimal step of the tracked particles, (b) in- 
clude the treatment of the longitudinal modes of oscillation at 0[(R/ L) 2 ], 
(c) evaluate the effects due to the beam structure. This study leads to a value 
of -0.04 for e in equation (1). 

> 16000 r 1 1 




AE in bar Me\ 

Fig. 5. Monte Carlo simulation for 510 MeV primary electrons impinging on the 
bar: AE per electron 



5 Measurements and comparison with expectations 

Several electron pulses were applied to the bar at various temperatures (264 K, 
71 K, 4.5 K) with DAQ collecting ADC data, thermometer data and number of 
electrons per pulse. The maximum amplitude of the first longitudinal mode of 
oscillation of the bar, B m , is related to experimental measured values according 
to the relation: 

B m = V max /{G\) 
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Table 2 

.Bo due to W=10~ 5 J and input values for the calculation in case of pure aluminum. 



T[K] 


m 


Am 


264 


0.96 


0.01 


71 


0.98 


0.03 


4.5 


1.16 


0.03 



Table 3 

Values of m fitting B m = m B t h- 

where V max is the maximum amplitude of the component of the signal at the 
frequency corresponding to the first longitudinal mode, obtained using Fast 
Fourier Transform procedures applied to the ADC stored data, G is the gain 
of the amplifier and A is the coupling factor defined in equation (3). 

Knowledge of ot(T) and c v (T) for aluminum is needed for computing the re- 
lation (2) at different temperatures: spline interpolations on data of ref. [19] 
(12 K< T <300 K) and the parametrization in ref. [20] (T <12 K) give a(T), 
while c v (T) is obtained by spline interpolations on values of c p reported in 
ref. [21]. Table 2 shows the calculated values of a, c v , B due to a deposition 
of 10 -3 J in a RAP-like bar made of pure aluminum. 

The correlations among measured B rn and estimated B t h (equation(l)) are 
shown in Fig. 6 and the values of the parameter m fitting the relation B m = 
m B t h are shown in Table 3. A systematic error of ~ 7%, obtained from the 
quadrature of the beam monitor (3%) and A determination (6%) accuracies, 
affects the measurements. 



6 Conclusions 

The results we obtained for the maximum amplitude of the fundamental mode 
of longitudinal vibration of the bar are in good agreement with the expecta- 
tions coming from the thermo-acoustic model describing the particle energy 
loss conversion into mechanical energy in the temperature range 270-4 K. This 
is the first time that experimental results on thermo-acoustic energy conver- 
sion in A15056 are obtained below 270 K using our technique. The result at 
T ~ 270 K improves the agreement with the model when compared to the 
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Fig. 6. Measured maximum amplitudes [lCr 1G m] of induced vibrations, B m (verti- 
cal), versus the expected ones, B t ^ (horizontal), for various beam intensities. Straight 
lines are the least square fits with equation B m = m B t h ■ 

findings of past experiments [3,22] based on the same technique. Amplitude 
measured values at T ~ 4 K are slightly higher than expectations: this is 
probably due to a different behavior of A15056 respect to pure aluminum or 
to an imperfect knowledge of the input parameters for the calculations in a 
region of temperatures where the ratio a/c v sensibly depends on T. 

We are greatly indebted to Dr. S. Bertolucci for the continuos support given 
to the experiment, to Drs. D. Babusci and G. Giordano for having derived the 
formal expression of the divergences used in the simulations of Appendix A 
and to Prof. M. Bassan for many useful discussions. We would like to thank 
Messrs. G. Ceccarelli, R. Ceccarelli, M. De Giorgi for the help given during the 
cryogenic operations of the detector and Messrs. M. Iannarelli and R. Lenci 
which helped with the experiment readout setup. 



A Monte Carlo simulation 



The energy in the normal modes of order n of a thin cylindrical bar (radius 
R much smaller than the length L) of mass M excited by a ionizing particle, 
can be calculated using the thermo-acoustic model [9] in terms of the sound 
source S, of the thermo- mechanical parameters of the bar and a geometrical 




T=71 K 
m = 0.98 
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factor G 



n ■ 



1 f L ■ G n \ 2 1 f L • G n dE\ 



where 7 is the Griineisen parameter, v s is the speed of sound, dE/dx is the 
specific energy loss. The geometrical factor G n , taking into account the path A 
inside the bar and the coupling to the eigen-modes of frequency u n , is defined 
as the path integral of the divergence of the displacement vector (normalized 
to the volume) u n : 

A 



In other terms, the energy for the n-th mode is obtained by squaring the 
amplitude of the vibration for a path of a single ionizing particle inside the 
bar: 

E - —^.A 2 

n ~2Mul n 



The corresponding expression when the bar is crossed by N particles is given 
by summing (incoherently) over all the amplitudes, 

N 

An = X] A$ ' 

i=l 



so that Eq. (A.l) can be written as: 

£ »^S(S^) 2 (A - 3) 



In the estimate of the expected amplitude when a beam of particles crosses the 
bar, two extreme approaches are possible: the specific energy loss is assumed 
to be constant all along the trajectory inside the resonant bar, so that it can 
be computed separately from G n ; or the real energy loss is evaluated for each 
infinitesimal step and is multiplied by the correct geometrical factor before 
summing all the amplitudes: 

A n =j:f(v-u n ) §di 

i=l{ aL i 
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The latter method requires a simulation of the elementary processes of en- 
ergy loss inside the material. The GEANT routines are able to track all the 
secondary particles produced in the interaction with the material of the gener- 
ated (primary) particle, and all the different physical processes occuring when 
the particles cross the matter are simulated, taking into account the appro- 
priate cross- sect ions for the different materials in the set-up, and the current 
energy and momentum of the particle: ionization, Coulomb scattering, and 
Bremsstrahlung for charged particles; pair production, Compton scattering, 
and photoelectric effect for photons. 

All the secondary particles produced in the interactions with the matter are 
followed until their energy falls below a given threshold. The size of the step of 
the simulation for each particle is automatically adjusted, in order to optimize 
the accuracy vs. the computing time, taking into account the set-up geometry 
and materials together with the kinematical quantities of the currently tracked 
particle. 

Such a simulation code is well suited to include the calculation of the geo- 
metrical integral in the tracking routines, together with the evaluation of the 
specific energy loss for the different physical processes happening inside the 
materials, at each step of the simulation. Moreover, since the contribution of 
each single particle in the beam has to be added in an incoherent sum, it is 
easy to extend the sum over all the secondary particles of the electromagnetic 
shower that develops inside the bar for each incident electron. 

The total excitation amplitude for the i-th primary electron in the beam is 
then evaluated adding the contribution of all the Ni secondary particles (j = 
is the contribution of the primary electron itself), at each of the Mj tracking 
steps in the simulation code: 

4 l) = EE 

j=0 k=l 

By summing together the contribution of all the steps of all the tracked par- 
ticles, quantities such as the total path length: 

j=0 fc=l 

and the total energy loss: 

j=0 k=l 
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can be calculated for each incoming primary electron. 

The integral in the expression of G n is numerically calculated for each step, 
assuming that the divergence of the displacement vector is constant (this is a 
good approximation for steps of ~ 1 mm): 

= J (V • u n ) —dl 3 = £ (V • u n (x jk )) ■ AE jk . (A.4) 

a . 3 k=l 



In order to compute the divergence V • u n (xjk), the exact solution of the 
normal modes of a finite cylindrical bar have been replaced by the approximate 
solution for an infinite cylinder due to Pochhammer and Chree (PC). This 
approximation holds for a thin cylinder, i.e. as long as a n = niij; < 1. The 
PC solution has been expanded to the order O[o%\ leading to: 



V • u n {x jk ) = 



V2a n 
R 



1 + 



a(l-a) 



2a 



a ( r 
1 + 2 [ an R 



sin [®n^ 



where a (=0.347 for isotropic aluminum at room temperature) is the Poisson 
module and r, z are referred to a cylindrical coordinate system having the z- 
axis coincident with the bar axis and the ^-origin located at one of the bar 
end-surfaces. 

In order to simulate a number of beam pulses on the RAP bar, the vibrational 
amplitudes in the different normal modes have been computed simulating a 
given number of electrons in the beam, with the nominal energy of 510 MeV 
and with the proper beam spread and divergence. This gives an estimate of 
the bar response for a given beam intensity. The following results refer to 100 
pulses of 10 electrons/pulse, hitting the bar at the center. The distributions 
of the A n amplitudes for the first four modes are shown in Fig. A.l. The odd 
modes have a non-zero average, while the the average for the even modes is 
consistent with zero. This is what we expect from the properties of the eigen- 
functions, that are symmetric (anti-symmetric) with respect to z = L/2 for 
odd (even) modes. Since the beam enters the bar at z — L/2, A ra=13j will 
have the same sign; while A n=2 ^... will change sign and have a null expectation 
value. 

In order to obtain an estimate of e (equation (1)) from this full simulation 
we compare the average value of Ai=-0.508 GeV/m, derived from the cor- 
responding histogram of Fig. A.l, against the value of A\ when computed 
Ola 1 ]: 

A 1 = ^(2cr — 1) < AE >= -0.530 GeV/m , 

Lj 
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Fig. A.l. Monte Carlo results for the first four modes of the amplitude: on the left 
the odd modes (n = 1,3,), on the right the even modes (n = 2,4) distributions (on 
top, bottom respectively). In the histograms all the events from 100 simulated shots 
of 10 4 electrons with E = 510 MeV have been accumulated. 

which is related to equation (2) in the hypothesis that all the energy lost by 
an electron (cfr. equation (4)) is deposited at the bar center. Thus we obtain 
e ~ -0.04. 
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